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Abstract 

We analyze the structure of running coupling corrections to the gluon production 
cross section in the projectile-nucleus collisions calculated in the Color Glass Condensate 
(CGC) framework. We argue that for the gluon production cross section (and for gluon 
transverse momentum spectra and multiplicity) the inclusion of running coupling correc- 
tions brings in collinear singularities due to final state splittings completely unaffected by 
CGC resummations. Hence, despite the saturation/CGC dynamics, the gluon production 
cross section is not infrared-safe. As usual, regularizing the singularities requires an in- 
frared cutoff AcoU that defines a resolution scale for gluons. We specifically show that the 
cutoff enters the gluon production cross section in the argument of the strong coupling 
constant (^^(A^qjj). We argue that for liadron production calculations one should be able 
to absorb the collinear divergence into a fragmentation function. The singular collinear 
terms in the gluon production cross section are shown not to contribute to the energy 
density e of the produced matter, which is indeed an infrared-finite quantity. 



1 Introduction 

In the recent years there has been a lot of progress in understanding running couphng cor- 
rections for small-x evolution equations. Running coupling corrections to the Balitsky-Fadin- 
Kuraev-Lipatov (BFKL) [1,2], the Balitsky-Kovchegov (BK) [3-7] and the Jalilian-Marian- 
lancu-McLerran-Weigert-Leonidov-Kovner (JIMWLK) [8-15] evolution equations were cal- 
culated in [16-19]. Knowledge of the running coupling corrections to the BK and JIMWLK 
equations significantly improved our ability to make precise predictions for the total cross sec- 
tions in deep inelastic scattering (DIS) based on the saturation/Color Glass Condensate (CGC) 
physics [3-15,20-30]. 

However, to improve the CGC predictions for other observables, such as the inclusive gluon 
or quark production cross sections, it is important to understand how the running coupling 
corrections enter the expressions derived previously for such cross sections [31-47]. It was 
originally proposed in [20,48] that the gluon production cross section can be described by the 
/cT-factorization formula 



where the gluon is produced with the transverse momentum p {pt = \p\), rapidity y in a 
collision with the total rapidity interval Y. (pp and (pt are the unintegrated gluon distribution 
functions in the projectile and in the target correspondingly, is the strong coupling constant 



/cy-factorization ([T]) has been known to work outside the saturation region [31-34, 49] (e.g. 
for proton-proton collisions pp). In [38,41,50] it was shown that the /cy-factorization formula 
([1]) holds for gluon production in proton-nucleus (pA) collisions and in DIS both in the quasi- 
classical limit of multiple rescatterings [22-25] and after including quantum BK/ JIMWLK 
evolution. This means that Eq. ([1]) works in the kinematic region inside the saturation region 
for the target and outside the saturation region for the projectile. (A proton/projectile can 
be defined as a nucleus with a much smaller saturation scale than in the target nucleus.) It 
is at present not clear whether Eq. ([1]) correctly describes gluon production in nucleus-nucleus 
collisions, i.e., in the saturation region for both the target and the projectile. 

In [38] the following relation was also derived between the forward amplitude of a dipole of 
transverse size x scattering on a nucleus at impact parameter b and rapidity y, N{x, b, y), and 
the unintegrated gluon distribution function 0(q, y) of the same nucleus: 



It would be very interesting to see if Eq. ([T]) would still hold when running coupling correc- 
tions are included. It would also be interesting to see whether the relation ([2]) would survive 
the inclusion of running coupling corrections. This would verify the validity of /cy-factorization 
formula ([1]) and the universality of the dipole amplitude as the correct degree of freedom for 
both the total DIS cross section and for the single gluon production at small-x. 

To answer the above questions one has to first understand how to include running coupling 
corrections into the gluon production cross section. This is a formidable task that we are not 
going to complete here. Instead we will be content with a first step in that direction by resolving 




(1) 



and Cf = {N^ - l)/2Nc. 
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some conceptual questions and by setting up a formalism for future calculations. Our main task 
is to identify the cancellation mechanisms that ensure the IR safety of the total cross section, 
as it governs all left over phase space cancellations for less inclusive observables and allows us 
to pinpoint the origin of eventual IR problems. 

Below we will analyze the running coupling corrections to the one-gluon inclusive production 
cross section. Our strategy is standard: following Brodsky, Lepage and Mackenzie (BLM) [51] 
we will begin by dressing all the relevant gluon lines and vertices by quark bubble chains, 
resumming powers of Nj. We will then replace 

Nf^-QTiP2 (3) 
to complete it to the full one-loop QCD beta-function 

ft ^ (4) 

Resumming all factors of (32 should replace all factors of the bare coupling constant in 
the original lowest order expression ([1]) by the physical couplingsll] 

The paper is structured as follows. In Sect. [2] we will set up the formalism for including 
quark bubbles in the final state (the state after the projectile's interaction with the target). 
We will then analyze the running coupling corrections to gluon production in Sect. [3l 

The main conceptual problem in including running coupling corrections for the gluon pro- 
duction cross section is with insertion of quark bubbles onto the outgoing gluon line (see Fig. [TT] 
below). Indeed, if one uses dimensional regularization and the quarks are massless all such bub- 
ble corrections are zero, since the on-mass-shell gluon provides no invariant momentum scale. 
This zero can be thought about as a "cancellation" of the infrared and ultraviolet divergences 
in the bubbles. The zero result is unsatisfactory, since it would leave the gluon production cross 
section proportional to a factor a^, the bare coupling, and thus unrenormalized. 

To resolve this problem we notice that collinear divergences, resulting from gluon splitting 
into two massless quarks (or gluons) come in at the same order in perturbation theory as 
running coupling corrections on the outgoing gluon line. To keep track of such divergences we 
regularize them with the infrared cutoff Acou. With the regulator inserted, we show in Sect. 
|3]that the bubbles on the outgoing gluon line generate a factor of running coupling as^^^oii) 



"'^We refer the reader to [52-55] and [56-58] for a discussion of the context and validity the BLM resummation 
procedure or its refinements, for example via the dispersive method. The first issue here is whether the quark 
bubbles faithfully trace all running coupling corrections and thus justify ([3]): The answer given in the references 
is that this is true in the "single dressed gluon approximation," i.e. if the process only involves one gluon 
line that needs to be dressed. The second issue is the all orders resummation involved in renormalizing the 
coupling: This is deeply connected with the large order behavior of perturbation theory and the notion of 
renormalons [59,60]. Despite the natural quantitative uncertainties inherent in trying to answer the question of 
what is the scale in the coupling, the BLM method should at least yield a qualitatively correct answer and most 
importantly, be able to tell perturbative from non-perturbative contributions. If one takes the most conservative 
and pessimistic view, one can always reexpand the resummed coupling back to its lowest scale carrying order 
and estimate the uncertainty induced by varying the renormalization scale "within reasonable bounds." Our 
expressions must match these lowest order corrections as exemplified for the running coupling corrections to the 
non-linear small- a; evolution in the calculation of the gluonic component of the NLO BFKL, BK, and JIMWLK 
kernels done in [61]. In [61] the result of a full NLO calculation confirmed the scale-setting for the running 
coupling which had been done in [16, 17] using the BLM prescription. 



2 



with the infrared (IR) cutoff in the argument. Therefore, the gluon production cross section 
has to come in with a factor of the couphng running with some arbitrary IR cutoff in the 
argument. Indeed the cutoff has a meaning of the resolution scale for gluons: once such a 
cutoff is introduced, coUinear qq pairs (and collinear gluon pairs) are also identified as gluons, 
generating extra diagrams contributing running coupling corrections (see Fig. [T3l) . which in the 
end renormalize gluon production cross section. We conclude that the gluon production cross 
section depends on the IR cutoff/resolution scale and is thus not infrared- safe, in agreement 
with conventional wisdom. 

Our argument shows that the result of the full calculation of the running coupling corrections 
would necessarily generate a power of ^^(A^qii). We expect that, again in line with conventional 
wisdom, such non-perturbative factors of the coupling could be absorbed into fragmentation 
functions if hadron production is the goal of a calculation at hand. The separation of the 
diagrams into the perturbative cross sections and the fragmentation functions could be done 
using the conventional approach (see e.g. [62,63] and references therein). It would result in a 
factorization scale dependence for both the perturbative cross section and the fragmentation 
function. 

Indeed one might worry that when Eq. ([1]) is applied to pA collisions, the unintegrated 
gluon distribution in the projectile proton (pp would depend on some non-perturbative scale 
A characterizing the proton, thus making the gluon production cross section in Eq. ([T]) non- 
perturbative even without the inclusion of running coupling corrections. This is a legitimate 
concern, though, for instance, it does not apply to the case of DIS, where the scale characterizing 
the projectile is the photon's virtuality ^qcd' which makes both 0p and the cross section 

in Eq. ([T]) perturbative before the inclusion of running coupling corrections. However, even in 
the case of DIS, running coupling corrections would still lead to the collinear singularities that 
we discuss below, making the gluon production cross section not infrared- safe. The problem of 
applying Eq. ([1]) to pA collisions is expected to be remedied in nucleus-nucleus (^^4) collisions 
or in pA collisions at very high energies, where the IR-safety of the gluon production cross 
section at fixed coupling is expected to be imposed by the saturation effects in both the target 
and the projectile (see [31-34]). While we can not prove that here, we also expect the collinear 
singularities discussed below to arise if one tries to calculate the running coupling corrections 
to the gluon production in AA collisions. 

We continue our discussion by examining an observable built from quark and gluon pro- 
duction cross sections which is infrared-safe in Sect. HI The goal is to show explicitly how 
the infrared contributions cancel to render the quantity infrared- safe. The observable is the 
energy density of the matter produced in a collision. In [64] the energy density of the medium 
produced in a heavy ion collision at the space-time point described by the proper time r, space- 
time rapidity r] and transverse coordinate b = (fci,fc2) was shown to be given by the following 
formula: 

^ M _ ^ [ d^pp, ( "^^^ I "^^^ I ] (5) 

' ' T J \d?pdrid?h d^pdrfd^h d'^pdrjd'^h \ 

Here ^2p^'^2{, is the number of gluons produced with transverse momentum p, momentum-space 
rapidity y and at transverse coordinate h. The quark and anti-quark production contributions 
(ppdridH ^"^^ d'^pdrid'^b added to the gluon one on the right hand side of Eq. ([5]). Note that 
the momentum space rapidity y gets replaced by the space-time rapidity r] in Eq. ([5]). The 
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above formula ([5]) is valid for sufficiently late proper times, r 1/ {pt) with {pt) the typical 
transverse momentum of the produced quarks and gluons. The corrections to it are suppressed 
by extra powers of r [64]. Eq. ([S]), along with our discussion below, is applicable to the particles 
produced at central rapidity, i.e., particles with small-x both in the target and the projectile 
wave functions. Eq. ([5]) is valid not only for heavy ion collisions, but for pp, pA and DIS 
collisions as well. 

In Sect. Hlwe show that the energy density ([5]) is indeed independent of the IR cutoff Acou 
by showing that terms in the gluon (and quark) multiplicity which depend of this cutoff and 
which are singular in the Acou — > limit actually vanish when inserted in Eq. ([5]) using the 
abstract structure of the equation. In Sect. [5] we demonstrate how these terms vanish in an 
explicit toy model: we show that coUinearly divergent terms have such a p^-- dependence, that, 
when inserted in Eq. (jS]), they integrate out to zero. Thus energy density is an infrared-finite 
quantity. 

We also argue that the running coupling corrections are included in the expression for the 
energy density ([5]) in a manner similar to the total cross section calculation [16-19], since both 
gluons and quarks are allowed in the final state for either the total cross section or for the 
energy density. 

As an aside we discuss the renormalization of the interaction of the projectile wave function 
with the target in Appendix |X] using a quark- ant iquark dipole as the projectile. We find how 
the running coupling corrections enter the Glauber- Mueller formula (lAip [65]. We obtain that 
one of the two couplings in the exponent of Eq. (lAip runs with the inverse transverse size of the 
dipole, while the other coupling runs with some non-perturbative scale characterizing nucleons 
in the nucleus, as shown in Eq. ( 1A8I) . 

We conclude in Sect. [6] by discussing our main results and potential future developments. 

2 Total Cross-Section: Cancellation of Fermion Bubbles 
in the Final State 

We begin by analyzing the running coupling corrections to the total cross section arising in the 
final state. Our discussion will apply to both the cross sections calculated in the quasi-classical 
approximation of McLerran-Venugopalan (MV) model [22-24] and to the cross sections where 
the corrections due to non-linear small-x evolution equations [3-15] are included. Just like in 
the calculation of the total dipole-nucleus cross section, the corrections due to small-x evolution 
come in through a sequence of longitudinally soft gluon emissions, as only gluonic corrections 
generate leading logarithms of Bjorken x (i.e. powers of as \nl/x) [1,2]. Our discussion below 
can equally be applied to the quasi-classical gluon production in the MV model, or to the gluon 
production enhanced by small-x evolution corrections, which appears as a quasi-classical gluon 
emission in one rung of the small-x evolution (see e.g. [30,38]). Common to both cases is an 
eikonal emission vertex for the gluons under consideration. 

In [16-18,66] the running coupling corrections to the non-linear small-x evolution equations 
[3-15,27] were calculated for the first time. The BK and JIMWLK evolution equations were 
written for the total cross section of a projectile (dipole) scattering on a target: therefore 
the optical theorem is applicable and allows to neglect all the final state interactions. This 
was used in [17, 18] to specifically discard the final state running coupling corrections. In [16] 
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Figure 1: The time-ordered picture of projectile-target scattering amplitude introducing nota- 
tions to be used below. All multiple rescatterings are denoted by a broad (blue) vertical line. 
The solid narrow vertical line denotes the final state cut. (Color on-line.) 

one-loop final state corrections were studied at order Nf and their cancellation was shown 
explicitly at that level. While we indeed expect also the all orders resummed final state fermion 
bubble insertions to cancel due to the optical theorem, we expect the explicit mechanism of this 
cancellation to be far from straightforward. Thus we set out to explore how the cancellation 
works explicitly. 

Before we begin, let us introduce the notations. We will consider a high energy proton- 
nucleus or DIS collision and will work in the center of mass frame. In this frame the incoming 
proton/gg dipole (henceforth referred to as the projectile) comes in from the light cone time 
of x+ = — oo moving along the The wave function of the projectile may contain extra 

gluonic fiuctuations. The projectile scatters on the target instantaneously at the light cone 
time x+ = 0. Similar to [38,40,67] we will denote the instantaneous interaction with the target 
by a broad (blue) vertical line, as shown in Fig. [1] The broad (blue) line comprises all possible 
multiple one- and two-gluon exchanges between the nucleons in the target and the projectile 
(similar to Fig. [15] in the Appendix These interactions of the target field with projectile 
constituents or the emitted gluon shown explicitly sum into path ordered exponentials along 
the projectile light cone direction. They are indicated for the emitted gluon (or any explicitly 
shown interacting quark in what is to follow) in form of an oval (red) mark. After scattering 
on the target the system evolves further towards the final state at light cone time x+ = -|-oo. 
This final state we will denote by a vertical solid line as shown in Fig. [1] and refer to it as the 
final state cut, as appropriate when talking about diagrammatic contributions to probabilities 
instead of amplitudes. 

The amplitude squared, using a similar notation, is shown in the center of Fig. [2l There 
the light cone time goes from x+ = — oo to x+ = +oo both in the amplitude and in its complex 
conjugate. In the case of the projectile dipole one has to sum over all possible emissions of the 
gluon by the quark and the anti-quark, both in the amplitude and in the complex conjugate 
amplitude. In case of the proton projectile the emissions off all of the valence quarks have to 
be resummed. To save space we will introduce an abbreviated notation, shown on the right 
hand side of Fig. [2l The disconnected gluon line implies sum over emissions of the gluon by 
all quarks in the projectile. The broad (blue) line indicates all multiple rescatterings in the 
amplitude. This notation also makes it easier to discuss our results in relation to any type of 
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Figure 2: On the left: the scattering amphtude squared for the diagram from Fig. [H On the 
right: an abbreviated notation for the same amphtude squared. 



(small) projectile. 

In the notation of the graph on the right of Fig. [2] the optical theorem can be formulated 
as follows: to calculate the total scattering cross section the dynamics between x+ = and 
x+ = -l-cxo can be ignored. For instance, if we wanted to find contribution of extra softer gluon 
emissions to the total cross section, only emissions at light cone time — oo < x+ < should be 
included both in the amplitude and in the complex conjugate amplitude. This indeed has been 
verified in the literature [38,68]. The same property applies to running coupling corrections. 
The quark bubbles dressing the gluon line at times < x+ < +00 do not contribute to the total 
cross section according to the optical theorem. Such corrections were indeed not included in 
the calculations of [16-18,66]. However, these corrections should be included if one is interested 
in gluon production cross section. To learn how to best organize a calculation that includes 
them, we will first see how they cancel in the total cross section. 

To set the stage, let us briefly recall how the cancellations for the cross section mandated by 
the optical theorem work at leading order. The diagrams to consider for the amplitude include 
all contributions from order through g^: 



(6) 



The optical theorem guarantees that at order as, the only diagrams contributing to the absolute 
value squared of this amplitude are 
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(7) 



This quite drastic reduction in the number of diagrams one needs to calculate to obtain the 
total cross section is due to the following set of separate cancellations: 
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That the cancellation pattern is separated out into three groups is dictated by the interaction 
of the produced gluon with the target as marked by the (red) dot on the interaction area: only 
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diagrams with identical interaction patterns can possibly cancel against each other. This will 
be the first guiding principle to identify potential cancellations below. The actual cancellations 
require further analysis of the integrals associated with the diagrams. If one restricts the 
sum over the final state phase space, the cancellations are no longer complete. Still, partial 
cancellations from the phase space regions that have not been excluded from the final state 
sum remain in effect and only become manifest if one groups the diagrams as shown in ([8]). 

We will not analyze all contributions exhaustively but rather focus on a number of subsets of 
the one and two loop fermion bubble insertions that cancel separately in the total cross section 
to illustrate and identify the basic mechanisms at work. This will lay the groundwork to allow 
us to precisely pinpoint the origin of IR divergences, understand how they cancel in the total 
cross section. It also allows to understand, for example, the cancellations of UV logarithms 
of the form ln(s//i) present in individual diagrams (with /i the renormalization scale, a UV 
quantity) from the final result. All these relationships are a prerequisite for a full calculation 
of the running coupling corrections that remains beyond the scope of the present paper. 

We begin with a single fermion bubble. The relevant diagrams with final state interactions 
(splittings and mergers) are shown in Figs. [3] and |H Due to light-cone time-ordered nature of 




Figure 3: One- loop quark bubble final-state corrections to the contribution of the diagram in 
Fig. [2] to the total cross section. Here the quark bubble does not interact with the target. The 
dotted vertical line denotes an intermediate state labeled by a. 

these graphs, it is more convenient to use light-cone perturbation theory (LCPT) rules [69, 70] 
to calculate them. To understand the cancellations we will concentrate almost solely on the 
energy denominators: apart from complex conjugation on the right hand side of the cut, the 
numerators of all diagrams in Fig. [3] are the same. At the same time direct application of 
the rules of LCPT as stated in [69, 70] appears problematic due to singularities in some of 
the energy denominators. One of the singular denominators corresponds to the intermediate 
state a shown in graph A in Fig. [3l The light cone energy of the state a is identical to that 
of the final (outgoing) state, generating a 1/0-type singularity coming from the corresponding 
energy denominator. This problem has been discussed before in [68] for small-x evolution, and 
our discussion here will follow the calculations performed in [68]. To tackle the problem of 
singularities one has to calculate the energy denominators in the diagram in Fig. [3]A starting 
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from the usual Feynman perturbation theory. 

In going from Feynman perturbation theory to LCPT one performs integration over the mi- 
nus components of the internal momenta in the diagram. A Fourier transform of a denominator 
of a gluon or quark propagator becomes 

d k i —ik-T f d k dku- — » ( k+X—+ '°„ X-L — k-x] r ^, > ^, \ i \^ 

= y p^yj^^ ^ ^ M(.,)e(t,) + o(-x,)o(-t,)i 

(9) 

which effectively puts the particle on mass-shell assigning the light-cone energy (the minus 
component of the momentum) to be equal to 

- ^- (10) 
Thus a particle propagating over a light cone time interval x+ brings in a factor of 

The incoming particle with momentum fc) brings in a factor of 

e-iSfe^+ (12) 

while the outgoing particle with momentum fc) brings in 

(13) 

as can be inferred from Eq. ([9]). 

With the above rules in mind, the relevant contribution of the diagram A in Fig. [3] is 

^ - / / (14) 



where the overall minus sign comes from taking into account two powers of i coming from the 
quark-gluon vertices. The omitted (numerator) part of the diagram A is the same for all three 
graphs in Fig. O the gluon in all three cases can only have transverse polarizations. In Eq. 0141) 
we have introduced light-cone regulators (5 > 0, which make the integrals convergent at infinite 
light cone time. (Indeed, for the initial state emissions the regulators should come in with the 
opposite sign in the exponent.) The notations used in Eq. f|T^ are explained in Fig. [3]A.: the 
quark and anti-quark lines carry momenta (A;i+,fci) and (/c2+;^2) correspondingly, while the 
gluon carries momentum (g_|_, q) = (A;i+, fci) + {k2+, ^2)- The quark-gluon vertices are assigned 
light cone times xi+ and X2+, as shown in Fig. [3j\. Due to theta-functions in Eq. (Q one can 
see that only the ordering < xi+ < X2+ < +00 used in Eq. ( IT^ is allowed, as g+ > for the 
produced gluon. 

Performing the integrations in Eq. (IT^ and expanding the result in the powers of 6 yields 

A — ^ 1 1 1 _^ ^^^^ ^^^^ 

25 Eq — Ek^ — Ek2 2 {Eg — Ek^ — Ek^Y 
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Noticing that the diagram A' in Fig. [3] is equal to a complex conjugate of the diagram A we 
obtain 



A + A' = A + A* 



1 



{Eg - - Efcj 



o(5). 



(16) 



{A* denotes a complex conjugate of A.) The dangerous 1/6 singularity is canceled in the sum 
of the two diagrams. Calculating the diagram B in Fig. [3] in a similar manner (or directly by 
using LCPT rules from [69, 70] since there is no singular energy denominators in this graph) 
one gets 



B 



1 



{Eg - Ek^ - Ek^] 



o{6). 



Adding Eqs. f[T6|) and flTTI) and taking the 5 — > limit yields 

A + A' + B = 0. 



(17) 



We have thus shown the cancellation of the diagrams in Fig. [3] for the calculation of the total 
cross section. 

Diagrams C and D (and C and D') in Fig. [Hare in a separate class of diagrams where the 



D 




D' 



C 



Figure 4: One-loop quark bubble final-state corrections to the contribution of the diagram in 
Fig. [2] to the total cross section: here the quark loop interacts with the target. 

quark loop interacts with the target. They do not have singular energy denominators. Applying 
the LCPT rules [69,70] one obtains for the relevant parts of those graphs: 



C 



1 

Ek^ + Ek^ — Eg 



(19) 



and 



D 



Eg — -Efci — Ek2 



(20) 
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such that 



C + D = 0. 



(21) 



Similarly 



C + D' = 0. 



(22) 



This accomplishes the proof of the cancellation of one-fermion bubble corrections in the final 
state for the total cross section. We observe that cancellations happen within classes of diagrams 
that share the same type of interactions with the target: diagrams with a gluon interacting 
with the target cancel separately from the diagrams with the quark bubble interacting with the 
target. 

To see if the pattern persists at higher orders, let us analyze the diagrams with two fermion 
bubbles. The different classes with separate cancellation patterns are shown in figs. [5l [6] and [71 
We begin with diagrams that have only the gluon interacting with the target shown in Fig. [5l 
Their calculation proceeds along the lines outlined above. Note, that the diagram E in Fig. [5] 
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Figure 5: Two-loop quark bubble final-state corrections to the contribution of the diagram in 
Fig. [2] to the total cross section. Here the quark bubbles do not interact with the target. 

is a square of the absolute value of the diagram A in Fig. [3j however, in order to keep all the 
finite pieces of diagram E one has to expand the expression for the diagram A up to order 5 
before squaring it (i.e., one order higher than shown in Eq. ( lT5l) ). since the 1/5 singularities 
make order 5 corrections in A contribute to E. The final result yields 



E 



452 ^E,-Ek,-Ek,[ 



1 

+ 4 



{Eg - Ek:, 

1 



E^ 



A 5 {Eg- Ek, -Ek,y {Eg - Eu, - Ek. 



{Eg 



Eki 



Ek2) {Eq — -Efc:? — E, 



k4) 



{Eg 



Eki 



Ek'iY {Eq — E, 



ks 



{Eg-Eu,-Ek,f{Eq-E, 



ki 



E, 



ki) 



o{6). 



(23) 
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For the other diagrams in Fig. [5] one gets 



{Eq — — Ek^y^ {Ef,^ + Ek^ — Ek^ — E, 



k2, 



(24) 
(25) 



G 



1 



9 

+ - 



S"^ {Eg — -Efci — -E'fca) {Eq — Ek3 — EkJ 

1 3 



i Ek, + Eu, + 3^fc3 + 3^fc, - AE, 
; 5 {Eq — Eki — Ek^Y {Eq — Ek-^ — Ek^)"^ 



Ek,f 



+ 



O {Eq — Ek^ — Ek^) {Eq — Ek^ 

9 1 

8 {Eq — Ek^ — Ek^)^ {Eq — Ek^ — EkJ ' 



{Eg - Ek^ - Ek^y {Eq - Ek.^ - Eki) 



(26) 



G' = G*{h,k2 ^ h,h). 
Adding the contributions up yields 

E + F + F' + G + G' = 0. 



(27) 



(28) 



Similarly one can demonstrate the cancellation of the diagrams with the quark bubble 
interacting with the target on one side of the cut and the gluon interacting with the target on 
the other side of the cut, as shown in Fig. [6j 



Figure 6: Two-loop quark bubble final-state corrections to the contribution of the diagram in 
Fig. [2] to the total cross section: here only one of the bubbles interacts with the target. 



H + I + J + K = 0. 



(29) 



(The same also applies to the mirror reflection of those graphs with respect to the cut.) 

Finally, the diagrams with quark bubbles interacting with the target on both sides of the 
cut shown in Fig. [7] also cancel: 
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M' 



Figure 7: Two-loop quark bubble final-state corrections to the contribution of the diagram in 
Fig. [2] to the total cross section: both bubbles interact with the target. 

L + M + M' = 0. (30) 

We have thus demonstrated cancellation of fermion bubbles in the final state up to the two- 
loop level for the total cross section. We conclude that such cancellations happen at all orders. 
As was noted before, due to the optical theorem, no final state interactions enter the total cross 
section. Alternatively one could argue that for the total cross section it is not important what 
happens to the gluon after the interaction with the target: it may split into a qq pair, like in 
Fig. [3]B, or the qq pair may recombine back into a gluon, as in Fig. [3JA.. Due to probability 
conservation, the contributions of those two events should cancel, as the gluon remains a single 
gluon with probability one. This conclusion is true for any number of bubbles on the gluon line 
- all of them should cancel after the sum over all cuts. 

The above cancellations happen only if the momenta of the appropriate quark and gluon 
lines are equal in all canceling diagramqj. For instance, the momenta ki, ^2, q along with ki^, 
/c2+, and g+ have to be taken at the same values for all three diagrams A, B and A' in Fig. [3] for 
the cancellation of Eq. f|T8|) to take place. However, the allowed values of ki and k2 are different 
in the virtual diagrams A and A' and in the real diagram B in Fig. [31 In fact, as one can easily 
see, the diagrams A and A' in Fig. [3] give an ultraviolet (UV) divergence after the integration 
over the momentum in the loop. At the same time the UV divergence in the diagram [3]B is 
cut off by the center-of-mass energy s in the scattering problem. Therefore the cancellation in 
Fig. [3] is not complete: in fact, instead of the zero on the right-hand-side of Eq. (fTSj) . one should 
get a contribution proportional to ln(s//i) with the renormalization scale (a UV cutoff). One 
would get a paradoxical result: quark loop of Fig. [3] would bring in a leading logarithm of the 
center-of-mass energy s (i.e., a logarithm which comes in only with one power of the coupling 
in front). This is indeed impossible, since it is well known that the leading logarithms of the 
center-of-mass energy s are given by the gluonic contribution only [1,2]. 

The resolution of this apparent paradox is in the fact that, while the particular set of 
diagrams in Fig. [3] does generate an additional leading logarithm of the center-of-mass energy 
s, such logarithms get canceled if we sum all the relevant diagrams, a possibility that has 
been encountered already in [71] . Consider diagrams in Figs. H] and [S] which exhibit the same 
problem. For instance, diagram C in Fig. H] is UV divergent, while diagram D in Fig. H] is 

^We thank Jianwei Qiu for asking one of us a question which led us to the understanding of the subtle effects 
discussed below. 
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Figure 8: One-loop quark bubble initial-state corrections to the contribution of the diagram in 
Fig. [2] to the total cross section. 



cut off by s in the UV. Since the diagram in Fig. HP comes in with the overall minus sign 
compared to the diagram in Fig. [3]B, as one of the quark-gluon vertices is moved to early time, 
we conclude that diagrams C and D in Fig. H] give a contribution proportional to — ln(s/ /i). In 
the deep UV limit that we are studying the interactions of the quark bubble with the target 
are the same as the interactions of the gluon with the target: in the UV limit the quark and 
the anti-quark in the bubble would be very close to each other in transverse plane. Hence the 
remaining factors in the diagrams C and D in Fig. Hlare the same as those in the diagrams A, 
B and A' in Fig. [31 Therefore the logarithms ln(s//i) cancel between diagrams [3] A, B, A' and 
diagrams mC, D. Similarly one can show that the logarithms ln(s//i) cancel between diagrams 
MN, O, N' and diagrams H] C", D'. In the end, all the extra logarithms of s cancel, just as in 
the example of [71] for leading-order evolution. 

In fact, the above observation applies to the leading-order small-x evolution leading to the 
JIMWLK and BK equations [3-15]. There, in the forward amplitude, the diagrams with the 
gluon interacting with the target have a UV divergence, which should be regulated by the 
center-of-mass energy s. The diagrams with the virtual gluon have a divergence regulated 
simply by a cutoff. In the full evolution equation the divergences cancel between the real and 
virtual terms. To really see this, one has to gather all diagrams that become similar to the 
above in the UV to argue that logarithms of s cancel as well. One might be tempted to argue 
that the diagrams like the one shown in Fig. [3] do not contribute to the forward amplitude. 
While this is true for generic momenta in the bulk of phase space, it is in fact not the case for 
the deep UV limit, where the transverse momenta of the quarks in the loop are large, making 
the light-cone lifetime of the quark bubble very short, shorter than the time interval between 
the multiple rescatterings in the target. Therefore such diagrams in this UV limit have to be 
included in the forward amplitude. However, due to the cancellation patterns observed above 
we see that all such deep-UV contributions cancel exactly (according to Eqs. ffTSl) . fl2T]) . fl22|) . 
etc.), with the net effect that the sums of diagrams in Figs. [3l HI etc., do not contribute to the 
forward amplitude, despite the fact that the individual diagrams would do so. In reality one 
has to be aware of this subtlety and remember that the diagrams under consideration do not 
contribute only because of an intricate cancellation. 

We will now apply the machinery developed above to the calculation of running coupling 
corrections to the inclusive gluon production. 
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3 Inclusive Production: Summing Fermion Bubbles on 
the Outgoing Gluon Line 



Now let us try to include running coupling corrections to the inclusive gluon production in 
the scattering of a small projectile on a large target. A diagram contributing to the gluon 
production cross section is shown in Fig. [9] in the notation introduced in the previous Section. 



Figure 9: Gluon production diagram for projectile-nucleus scattering using the notation intro- 
duced above. 

The cross on the gluon line indicated that we tag on that gluon, keeping its momentum fixed. 

To include the running coupling correction we have to again "dress" all gluon lines and 
vertices with quark bubbles. To begin with let us point out that t-channel gluon exchanges 
with the target get renormalized separately from the rest of the diagram: each t-channel gluon 
with transverse momentum I comes in with a factor of the bare coupling constant and, after 
being dressed with fermion bubbles, it acquires the standard geometric series which can be 
absorbed in the denominator leading to the physical running coupling constant 



The details of this calculation are shown in Appendix |Al where we show how running 
coupling corrections enter into multiple gluon exchanges. There we consider dipole-nucleus 
scattering in the quasi-classical approximation. The results of the Appendix \^ while directly 
applicable to DIS, can also be applied to gluon-target interactions shown in Fig. O Using the 
crossing symmetry one can show that the gluon-target interaction on both sides of the cut in 
Fig. [9] is equivalent to a gluon dipole-target interaction [30, 38] with the gluon dipole made 
out of the gluon lines on both sides of the cut. At large Nc the gluon dipole is equivalent to 
two quark dipoles considered in Appendix |Al We expect the structure of the running coupling 
corrections found in Appendix |X] for a quark dipole to be preserved for the gluon dipole beyond 
large Nc limit. 

It is interesting to note that the vertices where the t-channel gluon connects to the nucleons 
in the target come with a non-perturbative coupling, whose scale is determined by the typical 
momentum scale in each nucleon, i.e., by the scale of order of Aqcd- The vertices where the t- 
channel gluons couple to the projectile or to a gluon in the projectile's wave function come with 
the coupling at the scale of the gluon's transverse momentum and are, therefore, perturbative. 

The important message for the discussion below is that the multiple rescatterings renor- 
malize separately from the rest of the diagram. The only remaining factor of bare coupling 
comes from the gluon emission vertex. In Fig. [9] the gluon is emitted before the interaction on 
both sides of the cut. Since summation over emissions off of all quarks is implied in Fig. [HI the 
emission vertices are not shown explicitly. The two vertices give a factor of a^: our goal is to 
see how it transforms into a physical running coupling constant. 





(31) 
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The gluon propagators at the hght cone time before the interaction, x+ < 0, can be dressed 
by quark bubbles both in the amphtude and in the complex conjugate amplitude, as shown 
in Fig. [TUJ Both for Feynman diagrams and in LCPT such quark bubble chains would give a 




Figure 10: Initial-state running coupling corrections to gluon production cross section, r and 
r' label the momenta of the gluon line on the left and on the right of the cut. 



geometric series each. The two chains of bubbles would contribute a factor 

1 



1 + /52 In ^ 1 + /52 In 



(32) 



where r and r' are the transverse momenta of the gluon before interaction with the target in 
the amplitude and in the complex conjugate amplitude, as shown in Fig. [TOl 

At the same time the vertices of the gluon interaction with the target could also have 
additive quark bubble (vertex) corrections. Adding the contributions of such bubbles gives a 
factor 



1 + In — 



Q 

1 + In 



(33) 



where Q and Q' are some momentum scales characterizing the vertices (which, in the actual 
calculations are likely to become combinations of momenta of several gluon lines involved). 

The combination of all bubble corrections in Fig. [10] and the bare coupling coming from the 
emission vertex is 



a, 



l + a^f32 In^ 




l + a^/32 InSj 


1 + /32 In ^ 




1 + In ^ 



a. 



(34) 



One can see that the bubble corrections of Fig. [TO] still fail to turn the emission vertex coupling 
constant into the physical running coupling constant. The only remaining bubble corrections 
which may complete to a physical coupling constant are the bubble corrections placed on 
the outgoing gluon line in the final state. 




'iiir(^)Tr(^)w^inf^^ 



Figure 11: Final-state running coupling corrections to gluon production cross section. The 
gluon carries momentum q. 

The quark bubble chain corrections on the tagged gluon line are shown in Fig. [TT] To sum 
them up we will use the result of Sect. [2] demonstrating that the sum of such corrections with 
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the ones with the qq pair in the final state cancels. Therefore, the sum of all the diagrams with 
the gluon in the final state, like those shown in Fig. [TT], is equal to the negative of the sum 
of all the diagrams with the qq pair in the final state, as illustrated in Fig. O The diagrams 
with the qq pair in the final state (on the right hand side of Fig. [T2l) have the advantage of not 
generating any singular energy denominators, like the ones discussed in Sect. |2j 



gluon 

cuts 





quark 
cuts 




Figure 12: Diagrammatic representation of the consequence of the final-state interaction can- 
cellations observed in Sect. El 

While we can calculate the diagrams on the right hand side of Fig. [12] using LCPT, it is 
more straightforward to do so using the standard Feynman diagram technique. In Fig. [13] we 
show the amplitude with the qq pair in the final state and with the gluon line dressed by quark 
bubbles. Our discussion below will be restricted to the case when the gluon is emitted in the 
projectile before the interaction, as shown in Fig. [13] The conclusions we will obtain will be 
straightforward to generalize to the case of gluon emission after the interaction. 




Figure 13: A diagram with a qq pair in the final state and with running coupling (quark 
bubble) corrections included on the gluon line. The negative square of this diagram gives the 
right-hand-side of Fig. | 



Let us consider massless quarks in the bubbles. In order to regulate coUinear divergences, a 
study of which is our primary goal here, we will regulate the integral over the relative momentum 
of the qq pair in the final state with some infrared cutoff Acou. 

Completing the factor of Nf coming from the fermion bubbles on the gluon line to the 
full beta-function by the replacement Nf —QTif32 we obtain the following expression for the 
amplitude depicted in Fig. [13] (for similar calculations see [17,40,72]) 



A = - Y, M%q,q, 



4 TT a,. 



A=±l 



In 



e^* ■ k (1 - 2a + \a) 5„y 



a[l — a) g+ 



(35) 



The notations used above are explained in Fig. [131 the outgoing quark and anti-quark have 
transverse momenta ki and ^2, while the gluon has 4- momentum g^. The longitudinal fraction 
of the gluon's momentum carried by the quark is 

fci+ _ 



a 



1+ 



2+ 



(36) 



16 



In Eq. ( |35|) we have defined 



k = {l-a)ki-ak2. (37) 

The quark and anti-quark helicities are labeled a and a' and the gluon polarization vector is 
= (0,0, e^) with = (l + zA)/y2 and A = ±1. In dia gram shown in Fig. [13] only transverse 
gluon polarizations contribute. The gluon carries color a and t"' is the color matrix: summation 
over repeated indices is implied. 

Remembering that outgoing quark and anti-quark are on mass shell and massless we write 

<l' = — (38) 
a a 

with a = 1 — a. 

The two-dimensional vector quantity M°-{q, q^) denotes the part of the amplitude including 
the gluon emission and its interaction with the target. The fact that M°-{q,q^) depends only 
on q and g+, and is independent of (or, equivalently, k), is most easily seen if this quantity 
is calculated in LCPT. There, all energy denominators in M"-{q,q+) would not include any 
fc-dependence since iVf (q,g+) includes only the initial state dynamics and k is the final state 
momentum. The factor of the coupling constant g due to gluon emission by the projectile is 
not included in M°-{q,q^) and is shown explicitly in Eq. fl35l) as one of the two factors of g 
giving in the numerator. 

As q"^ = k'^/aa > 0, and, since later on we will regulate the coUinear divergence at = 
with a cutoff Acoii, we drop the ie regulators in Eq. (l35i) obtaining 

A=y M'^iq, q^) . ^ , ^ ^ ' ^ {I - 2a + Xa) S.y ^ ^^^^ 



til ^ 1 + «m/52 ln(-±) 

L V 'MS y 



As the gluon is space-like, q^ > 0, the amplitude in Eq. fl39p has an imaginary part corre- 
sponding to the qq pair in any one of the fermion bubbles dressing up the gluon line going on 
mass shell. However, this imaginary part does not contain any ultraviolet divergence and hence 
does not contribute to running coupling corrections. We will neglect it along with the 5/3 in 
the denominator of Eq. fl39|l and rewrite the amphtude as 



A = 4n^ y M^{q,q,).e^^t^'-—^i^^^^^±Ml,^. (40) 

q-^ aaq+ 



A=±l 



Multiplying the amplitude fl40p by its complex conjugate and summing over colors, quark 
helicities a, a' and Nf quark flavors (assuming all flavors are massless) yields the amplitude 
squared 



\A\' = (4vr)^ y M\q,q^).e'M-{q,q^)-e''* 

A,A'=±1 \ Q / 



q\ 



X e^* -k e^' ■ [(1 - 2a)2 + A A']. (41) 

To calculate the contribution to the gluon production cross section of diagrams with the 
outgoing gluon line dressed with quark bubbles to all orders (but excluding the contribution 
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of the bare outgoing gluon line) we use the identity pictured in Fig. [12) we flip the sign of the 
amplitude squared in Eq. (1411) and integrate over k and a. The contribution to the total cross 
section of the square of the diagram in Fig. [T3] comes with the integration measure 

(fki dki+ (fk2 dk2+ d?q d?k dq+ da g+ , , 

2(27r)3 2(27r)3 " 2(27r)3 2(27r)3 ' ^ ' 

Defining the gluon's rapidity y by 

^ = lln^ = ln^ (43) 

with Qt = \q\ we write the contribution of the diagram on the left hand side of Fig. [12] to the 
multiplicity of produced gluons as 



d^qdy J 4(27r 



^6 



1 

Jdaql\A\\ (44) 







The integral over k in Eq. (I44p comes with an infrared regulator Acou: the domain of integration 
is restricted to k^ > Acon- On the other hand, when arriving at Eq. f[35l) one effectively integrates 
over momenta in the (uncut) quark bubbles up to the ultraviolet (UV) cutoff yUMS". Therefore, 
to keep the integration over kx in the cut bubble in Fig. [13] in agreement with the similar 
integrals in other (uncut) quark bubbles and thus to keep Fig. [T2] valid by restricting all quark 
bubbles to the same phase space, we will introduce this UV cutoff and put another restriction 
on the range of integration: kx < fijjig ■ It is important to stress once more that Eq. (jUj) does 
not contain the contribution of the bare outgoing gluon line. 

Indeed the equality in Fig.[T2]is true without any infrared cutoffs on k^p integration. In other 
words, the integral over k of the cut bubble on the right of Fig. [T2] is, in general, not restricted 
to kx > Acou. One can think of the right hand side of Fig. [T^] as consisting of the sum of the 
/ct > Acou and kx < Acou contributions (with an overall minus sign). In Eq. (jUj) we only keep 
the fcr > Acou contribution. Due to Fig. [12] it is equal to the sum of gluon production and the 
qq production with kx < Acou- Now the meaning of the multiplicity in Eq. (jH]) becomes clear: 
it is the net number of produced gluons and coUinear qq pairs with k^ < Acou- For simplicity 
we will refer to the obtained multiplicity distribution as "gluon" multiplicity below, but we will 
always keep in mind that coUinear qq pairs with kx < Acou are included in it as well. 

Replacing A''^ —611^2 in Eq- (1411) . plugging the result into Eq. (|44|) we integrate over the 
angles of k and using Eq. (j38|) for get 



M\q,q^)\'jda | ^4 (^^^^|^) ' ^ 4 [(1 - 2«)^ + 1]. (45) 

ALn 



d^qdy (2 7r) 

It is interesting to observe that the argument of the running coupling in Eq. fH5|) is similar 
to that of the coupling in the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) evolution 
equation [73-75] as was calculated in [76]. 
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Using 



■^QCD 

we integrate Eq. fHSl) over A;|,. The result of the integration yields 

1 



(Pqdy (2 7r)4 



bubbles _ - " ,^a^^^ ^^)|2 j _ 2^)2 ^ ^] |^ 

L 



coll \ _ ^ I ^MS 

etc 



(47) 



The part of the a-integral in Eq. (l47l) involving the factor of a a in the argument of the 
coupling constants contributes only to fixing the constant under the logarithm. In the spirit 
of BLM approach [51] we can expand the coupling constants in Eq. ( 1471) into geometric series, 
integrate each (of the first few) terms over a, and then resum back the series, a a terms in the 
argument of the couplings will only contribute constants to the terms in the series: since we 
are not keeping track of constant terms we will simply neglect them. Integrating the prefactor 
in Eq. fl47p over a and neglecting other irrelevant constants in the argument of the second 
coupling we write for the contribution to the produced gluon multiplicity 

HN^ 1 

= (2^ l^"'"' '^'1' - "J ■ 

Finally, adding the contribution of the bare outgoing gluon line 

= = p^|M"(,.,,)Pa„ (49) 

(with the bare coupling constant g) to the contribution of Eq. (148|) we obtain the following 
expression for the produced gluon (and collinear qq pairs) multiphcity with the running coupling 
corrections included: 

(Note that Eq. fl50p can be obtained directly from Eq. (HSj) by setting the upper limit of the k- 
integration to infinity. As follows from our derivation above, doing this would modify Eq. psl) 
to include the contribution of the bare gluon line. Similar observation has been made before 
for other production processes in [56,77].) 

From Eq. ( l50l) we can draw the following conclusions: 

• First of all, we have shown that the bare coupling constant stemming from the ver- 
tices of gluon emission by the projectile in the amplitude and in the complex conjugate 
amplitude gets renormalized into a physical coupling constant. We have thus outlined 
the procedure of how running coupling corrections should be included into the gluon 
production cross section. 
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• Unfortunately the scale with which the coupling runs in Eq. (!50|) is some arbitrary infrared 
cutoff Acou which we introduced to regulate collinear divergences. This makes the gluon 
production cross section (or, equivalently, gluon multiplicity) either perturbative or non- 
perturbative depending on the choice of the cutoff Acou corresponding to the resolution 
scale for gluons and qq pairs. Strictly speaking this means that the gluon multiplicity is 
not infrared-safe, in accordance with conventional wisdom. 

Indeed the dependence of the gluon multiplicity on Acou is already weakened by that fact 
that it is absorbed in the argument of a running coupling constant. The Acou-dependence could 
be further reduced by performing a Sudakov-like resummation [78] of gluon splittings. If the 
goal of a calculation is to find the cross section for production of hadrons, it may be possible 
to absorb the result of such resummation into a fragmentation function. However, performing 
a Sudakov resummation on top of the calculation of the running coupling corrections appears 
to be rather involved and is beyond the scope of this work. 

The resolution scale Acou generates an IR cutoff on kx of the produced quarks in Eq. (H5l) . 
Due to Eq. fl38|) this translates into a cutoff on gluon virtuality q^ , which is equal to the invariant 
mass of the qq pair. Neglecting a a we could say that qq pairs with the invariant mass q^ < A^^^ 
are included into our definition of gluon multiplicity distribution. 

However, one should be aware that our calculation is reliable only as long as A^^^ remains in 
the perturbative domain and that, if one is forced to take A^^^ closer to Aq^j-, on phenomenolog- 
ical grounds, the only consistent way to take into account the non-perturbative contributions 
- and the nontrivial true dependence on A^^jj in the non-perturbative region - is to introduce 
a factorization procedure that subsumes these contributions in a fragmentation function whose 
^cou dependence can only be inferred with non-perturbative means or by comparison with ex- 
periment. This introduces a factorization scale /i^ between perturbative contributions above 
it and non-perturbative ones below. The former are reliably reproduced by our calculation if 
we substitute the perturbative factorization scale jj? in place of A^^jj in the abovejfl For heavy 
quark-anti-quark bound state production the factorization scale is chosen to be equal to the 
lowest possible value of the invariant mass of the qq pair [79], which is not too different from 
the cutoff on the invariant mass obtained above. Still, the introduction of fragmentation func- 
tions is necessary for consistency even in such a case where the non-perturbative contributions 
become (almost) trivial. 

In practical phenomenological calculations of hadron production one often chooses the fac- 
torization scale fi ^ pt with pt the transverse momentum of the produced hadron. In this 
case, Eq. ( l50i) should be convoluted with the fragmentation functions and Acou in it should be 
replaced by the factorization scale fi. If one puts ^ ^ pt one would then get as(pf.) in Eq. (l50l) . 
While in an all-order calculation the hadron production cross section should not depend on the 
choice of the factorization scale /i, in practical applications putting /i ~ pj- in our result would 
lead to some additional p^-dependence. 

■^In such a setting, one has to keep in mind that the factorization scale, while infrared from the perspective 
of the perturbative contribution, is in the UV from the perspective of the modes entering the operator defining 
the fragmentation functions [62, 63] and may be thought of as its renormalization scale /i, which, in general, 
is different from full UV renormalization scale fij^ for the complete graph. Fragmentation functions depend 
both on the factorization scale n and on the collinear IR cutoff AcoU- The perturbative part of the cross section 
depends on the renormalization scale /if^ and on the factorization scale fi. In the full hadron production cross 
section the /i-dcpcndcncc cancels so that it only depends on fJ-jj^ and Acoii- 
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4 Collinear Singularities and the Energy Density 



Our goal now is to verify that while gluon (and quark) multiplicity is not infrared- finite, the 
energy density of the medium produced in the collision is infrared-finite. The energy density e 
is given by Eq. ([5]). With the gluon multiplicity being infinite, there are at least two ways the 
net energy density could stay finite: 

(i) it may be that the sum of gluon and quark multiplicities (net parton multiplicity) is finite, 
or 

(ii) it may be that the divergence in the sum of gluon and quark multiplicities vanishes after 
the integration over the transverse momentum in Eq. (JSj). 

To calculate the energy density using Eq. ([5]) we need the expressions for gluon and quark 
multiplicities. For the gluon multiplicity we have Eq. ( l50l) . To calculate the contribution of the 
diagram in Fig. [13] to the quark production cross section, we simply have to take the expression 
in Eq. psjl . fiip its sign back (to undo the sign fiip shown on the right hand side of Fig. [T2l) 
and instead of integrating over k and a we would integrate over momenta of the anti-quark 
keeping the transverse momentum and rapidity of the quark fixed. We get 



(Pp dy 



(2vr) 



d'^ki d'^k2 



dki+ dk2+ 



5{p — fci) (5 y — In 



ki+V2 



"IT 



\M'^iq,q+)\'f{k;a,a) 
(51) 



were we have defined 



f{k] a, a) = 6 7r/?2 



as{k^/a a) 



k^[{a-ay + 1]. 



(52) 



The produced quark's rapidity is given hy y = ln(p+ \/2/pt), same as Eq. (143!) for the gluon ra- 
pidity. Note that f{k; a, a) is symmetric under the interchange a ^ a: f{k; a, a) = f{k] a, a). 

Strictly speaking to calculate quark production we would have to replace —6 vr in Eq. (l52l) 
(and, consequently, in Eq. ( 1511) ) back with Nf. However we will leave (32 intact in Eq. (l52l) . 
By doing so we will indeed include the cut gluon bubble into Eq. (|5T1) making the expression 
in Eq. (|5T|) include both quark multiplicity and the multiplicity of gluons coming from a cut 
gluon bubble correction. Indeed if we want to see cancellation or disappearance of collinear 
singularities, we need to sum both (anti-)quark and gluon contributions. Hence, in our notation, 
the "quark multiplicity" would include both quark and gluon multiplicity coming from the cut 
bubbles. 

Eliminating the delta-functions in Eq. (ISTl) we write 



dm 



d'^k2 I da 



f{ap — a k2] a, a) 



(53) 



were, in accordance with Eq. (I36l) . we have used 



a 



P+ 

p+ + k2+ ' 



(54) 
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The anti-quark multiphcity is obtained from the multiphcity of the quarks in (l53l) by re- 
placing and ^2 — > ki. We get 



d?pdy {211)' 



d^ki I da 



Pt 



a 



V2 



f{ap — a ki] a, a) 



(55) 



With the help of Eqs. (145|1 and (152|) we rewrite Eq. (l50l) as 



d^pdy (2 7r) 



d k da 



M" p 



Pt 



/(fc; a, a) 



(56) 



where now the integral over kr is not bounded by the cutoff yU^jg from above. As was discussed 
in the previous section this insures that the contribution of the bare outgoing gluon line is 
included in Eq. ( !56|) as well. 

Defining the net parton multiplicity as the sum of quark and gluon contributions 



dN^ 



dN'^ dm 
+ 



d^pdy d^pdy d^pdy d^pdy 
we write, using Eqs. ( 153!) . (!55|) and (!56ll 



(57) 
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f{ap — ak;a, a) 



(5^ 



where we have relabeled the integration variables ki and k2 in Eqs. ( I55ll and (1531) as fc. 

To analyze coUinear divergences in Eq. ( l58l) it is convenient to redefine transverse momenta 
in the arguments of the functions / in the second and the third term in the curly brackets in 
Eq. f l58|) . We replace ak — ap—>^k in the second term and ak — ap—>k in the third term 
to obtain 



^jyparto 



d'^pdy (2 7r)4 



d'^k / da f{k; a, a) 
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tt ' aV2 
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a2 



fp + k pTe^ 



a ' aV2 



(59) 



With the help of Eq. fl59|) we clearly see that the coUinear divergence at = does not get 
canceled between the terms in the curly brackets. 

One might imagine some special "fine tuning" choices of amplitudes M (e.g. 



Pt 



i/p 
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without any rapidity dependence) for which such cancellation would be possible: however, even 

2 



M 



V2 



1/p: 



and cancellation 



at the lowest order in the coupling one has [31,34] 

does not happen in Eq. ( |59l) . Higher order rescatterings or other corrections are not likely to 
modify this scaling of the amplitude squared for sufficiently high pt- Hence the special choice 



of 



M 



P 



V2 



~ 1/Pr Pt appears to be not achievable in QCD. 

Therefore, the possibility (i) outlined above is ruled out: the net parton multiplicity still 
contains coUinear divergences. 

Before we proceed it is important to point out that, while similar to Eq. (145|1 the fc-integrals 
in the second and third (quark and anti-quark) terms in the curly brackets in Eq. (l59l) are cut 
off by in the ultraviolet, the A;-integral in the first (gluon) term in the curly brackets is not 
bounded in the UV. This is done to include the contribution of the bare outgoing gluon line, as 
was discussed in the previous Section and noted again in this Section after Eq. fl56l) . Indeed the 
UV behavior of the integrals is not essential for our conclusion of non-cancellation of coUinear 
divergences in Eq. (!59l) . but it will help us understand the physical picture later on. 

To verify the cancellation of coUinear divergences in the energy density, we first combine 
Eq. (15]) and Eq. ([STD to write 



e(r,?7, b) 



d ppt 



d'^p drj 



(60) 



In arriving at Eq. ( 160|) we have assumed the target (nucleus) to be cylindrical oriented along the 
collision axis. This allowed us to simplify the calculation by replacing d% in the denominator 
by the transverse cross-sectional area of the nucleus S±. This simplification does not affect 
coUinear divergences. 

Substituting Eq. (j59l) into Eq. ( !60l) we write 



e{T,ri,b) 



1 



TS±{2 7r 



^ I d^ppr I d^k I da f{k;a,a) 



M" p 



Pt e ' 



a a v2/ 



a2 



M 



„ ( p + k Pt e"^ 



(61) 



Again, the coUinear singularity is at fc = 0. Near the singularity we can put fc = in the 
arguments of in the second and the third terms in the curly brackets of Eq. (l6Ti) . The 

(potentially) singular part of the energy density is 



1 



^sing 



tSi (2tt)^ 



d^p Pt / d^k I da f{k;a,a 



M" p 



Pre'' 
V2 



a^ 



a aV^J 



a^ 



p Pt e ' 



a' ay/2 



(62) 



Rescaling pt by a in the second term in the curly brackets and by a in the third term in the 
curly brackets in Eq. fl62l) we see that the coUinear singularities cancel in the energy density! 
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Therefore, the option (ii) outhned above is correct! Energy density is independent of 
colhnear sphttings, while parton multiphcity is not. 

Let us point out that, as was observed before, the /cj^-integrals in the second and third 
terms in Eq. (I62l) have a UV cutoff /ij^, while the first term is not limited in the ultraviolet. 
Therefore, strictly speaking Eq. fl62|) gives 



vr 



^sing 



(fp Pt 



da f{k 



a, a) 



2 

MS 



Pre'' 
V2 



(63) 



(For simplicity again we do not distinguish between and, say, /iMs/'^O While the term 
in Eq. fl63p has no colhnear divergences with fc^ integral running over the UV modes, it is 
nevertheless non-zero. Integrating over kx and a we obtain 



^smg 



d^p Pt a. 



Pt e ' 



(64) 



This is indeed the contribution of the bare outgoing gluon line to the energy density. 

The physical picture behind the transition from Eq. fl62l) to Eq. flM|l is clear. When a gluon 
splits into a colhnear qq pair (or into a pair of gluons) the energy density deposited into a 
region of space would not change (in the strictly colhnear limit). Therefore, energy density 
is not affected by colhnear divergences/splittings. What remains of the energy density fl62|) . 
which included all real and virtual splitting, is the energy of the original gluon in Eq. flMl) . as 
it should be: none of the colhnear splittings would change the energy density. 

Of course one may be worried that the energy density in Eq. (l64l) contains a bare coupling 
constant instead of a physical coupling a^. This bare coupling gets renormalized in a manner 
similar to what was done in [16, 17] for the total scattering cross section: one has to add to the 
bare gluon term a diagram with a quark bubble (to be completed to the full beta-function by a 
similar gluon bubble along with other gluonic corrections), where the quark bubble both begins 
and ends in the initial state. This is shown in Fig. [HI As was demonstrated in [16,17] the quark 
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Figure 14: Addition of the diagram with the quark bubble (on the right) needed to renormalize 
the energy density e discussed in the text. The diagram on the right was also needed in [16,17] 
to take into account the running coupling corrections to the total cross section, as calculated 
using small-x JIMWLK and BK evolution equations. 



bubble diagram contains a UV divergence coming from the region of the transverse coordinate 
integral where the quark and the anti-quark are on top of each other. This divergence was 
instrumental in renormalizing the total cross section in [16,17], and is likely to renormalize the 
energy density in Eq. (IMl) . 
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5 Toy Model 



In this Section we will use a toy model for the amplitude squared to illustrate how the collinearly 
divergent expression for the net parton multiplicity (159|) leads to a finite energy density after 
the integration in Eq. (!60!l . Using the same approximations and variable redefinitions that led 
to Eq. ( !62l) we write the singular parts of the gluon and quark spectra as 



(Ppdy (2 7r) 



d k da 



MUp, 



(2vr) 



|M'^(p,P+)r«.(Acou), 
(65) 



d^p dy 



(2vr) 



d^k 



da — 
a^ 



a a\/1) 



f{k;a,a), 



(66) 



d^pdy 



1 



(2vr) 



d^k I da — 

a^ 



a aV2j 



f{k;a,a). 



(67) 



For simplicity let us write 

f{k; a, a) = Gn P2 



4 [{a-af + ll 



(68) 



which is obtained from Eq. ( l52i) by dropping a a in the argument of the coupling constant. 
Indeed, in all of the above calculations that factor was neglected in the end: we might as well 
neglect it right away here. 

Using Eq. (!68l) in Eqs. (l66l) and (1671) we perform kx integrals (remembering that Acou < 
/ct < /^Ms) obtaining 



d^p dy d^p dy (2 tt) 



da — [{a — a) + 1] 



P Pre 



a a\f2 



(69) 



Now, as a toy model, and in the spirit of saturation approach let us take 

1 



(2 7r)^ 



■jki ,pt < Qs 



(70) 



with some proportionality coefficient, which we put to one for simplicity. 

With the amplitude squared given by Eq. (170|) substituted into Eq. (!65|) we get 



d'^p dy 



a 



:(ALl) 



jr ,Pt > Qs 



(71) 
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for the number of produced gluons. 

Substituting Eq. fl7UI) into Eq. (1^9]) and integrating yields 



7 



1- ,pt> Qs 



d^pdy d^pdy L V coiv 2 , 3 3 1, i2p^ „ 

The term containing in Eq. (1721) is the contribution of the bare gluon, which is independent of 
the collinear cutoff Acoii. Therefore, since we are interested only in the cancellation mechanism 
for collinear divergences, we will drop this term in the further calculations. 

Using Eqs. (jTB and (j72l) in Eq. ( !57j) we obtain the Acou-dependent part of the net parton 
multiplicity 

,^partons [ ,PT > Qs 

^^"^ ^ «.(ALO (73) 



PtQ^ 2Q4 + q4 ill -h ,Pt<-^s 



d^pdy 

Now one can easily see that plugging Eq. (I73|) in Eq. (|60|) and integrating over p^ would give 
zero contribution to the energy density: 



The coUinearly divergent piece given in Eq. (1731) integrates out to zero! This is indeed due to 
the fact that dN^^^°^^ / d'^p dy is negative for some values of in the range of pr < Qs- 

What we have learned from the above toy model is that the collinearly divergent part of the 
net parton multiplicity dN^^'^^°'^'^ / d'^p dy , while non-zero, comes with a pr-dependent coefficient, 
which integrates out to zero when weighed with pt and integrated over p'^ to obtain the energy 
density using Eq. (!60|) . This clarifies the question of how exactly the energy density remains 
finite while dN'^^^^°'^^ / d^pdy is not. 



6 Outlook 

In this paper we have outlined the inclusion of running coupling corrections into the multi- 
plicity distribution of the produced gluons dN^ / d^pdy for the scattering of a small projectile 
on a larger target. We found that both the running coupling corrections and the collinear 
singularities enter at the same order in the coupling. Therefore, inevitably the gluon multi- 
plicity distribution turns out to be dependent on the infrared cutoff Acou which we inserted 
to regulate the collinear divergences. This result is indeed in agreement with the conventional 
understanding of gluon (or any other massless parton) production. However, when we resumed 
running coupling corrections (defined as powers of Nf with Nf later completed to the full 
bet a- function) to all orders, it turned out that the IR cutoff Acoii sets the scale of one of the 
factors of the running coupling constant, as shown in Eq. (I50p . This would indeed make gluon 
multiplicity distribution not infrared safe. For hadron production calculations one would expect 
that the contribution of the collinear singularities could be factored out into the fragmentation 
function. 
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We have then proceeded to analyze the effect of coUinear singularities on the energy density 
of the medium produced in the collision as defined by Eq. above (see [64]). In Sect. H] 
we show that the energy density e is independent of Acou- Therefore the energy density is 
infrared safe. This result is easy to interpret physically, as coUinear splittings, where a gluon 
splits into two gluons (or into a qq pair) both of which are flying in the same direction, are 
not going to change the amount of energy density deposited into a given region of space. We 
have also clarified how the parton multiplicity in Eq. ([5]) could be divergent due to coUinear 
singularities but give a finite energy density after pt weighing and integration in Eq. ([5]). It 
turns out, as was shown in Sect. [5l that the coUinear ly divergent part of the parton multiplicity 
comes in multiplied by a function of pt which, when weighed with p^ and integrated over all p"^, 
integrates out to zero. It is interesting to note that the energy density in Eq. ([5]) is proportional 
to another infrared-safe quantity called the energy flow, which was defined in [80] as a potential 
signal for the experimental detection of the BFKL evolution in deep inelastic scattering (DIS). 

CoUinear singularities modify the proof presented in [64] of the impossibility of perturbative 
isotropization (and thermalization) in heavy ion collisions which led to Eq. ([5]). To account 
for the coUinear singularities one has to separate each function fi{p'^,p ■ p',p'^,pt) employed 
in [64] (with k used in place of p in [64]) into coUinearly divergent and finite at = (or 
p'^ = 0) parts. For the singular part one has to introduce an IR regulator Acou as a mass in the 
denominator of the propagator of the outgoing gluon (or quark for the part of the argument 
made for the quarks). As was shown above in Sections H] and [5l the terms in e which are singular 
in the Acou — > limit will vanish after p'r-integration. For the remaining coUinearly-finite part 
of the function fi{p'^,p ■ p',p''^,Pt) the proof continues as was outlined in [64]. The conclusion 
of no perturbative isotropization is therefore not affected by coUinear singularities. 

While in this paper we have outlined the inclusion of running coupling corrections into the 
gluon production cross section, the full calculation still remains to be performed. It would 
indeed be important for improving our understanding of which part of the contributions can be 
expected to be made infrared safe by the presence of the saturation scale. As long as the other 
contributions can be treated with factorization methods, albeit at the price of having to intro- 
duce non-perturbative quantities like fragmentation functions, one arrives at a perturbatively 
consistent calculational framework. In such a framework, knowledge of the scales of the running 
couplings in the unintegrated gluon distribution functions which enter the fcr-factorization for- 
mula for gluon production in Eq. ([1]) [20,48], which was shown to be valid for gluon production 
in pA collisions and in DIS [35,38,41], is of paramount importance. It would allow one to make 
predictions for hadronic spectra and multiplicities in nuclear and hadronic collisions based on 
CGC physics with a much higher precision than was ever possible before. Such results would 
be extremely valuable for the upcoming LHC and EIC/eRHIC experimental programs. 
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A Running Coupling Corrections for Mueller- Glauber 
Multiple Rescatterings 

Here we will discuss running coupling corrections to multiple exchanges of t-channel gluons. We 
will consider scattering of a qq dipole on the nucleus in the quasi-classical approximation with 
two-gluon exchange interactions between the dipole and each nucleon, as was first calculated 
in [65]. This approximation is equivalent to the classical McLerran-Venugopalan model [22-27]. 
To concentrate on the running coupling corrections to the t-channel gluons let us take a quark- 
antiquark dipole as a projectile. The forward amplitude for multiple rescatterings resummed 
in [65] is shown in Fig. [151 As was demonstrated in [65] a sum of such multiple rescattering 




Figure 15: Glauber-Mueller multiple rescatterings of a quark-antiquark dipole on a nuclear 
target. 



diagrams gives the forward scattering amplitude 

N{x,b) = l-exp(-7ral^pTib)x'ln---]. (Al) 



The above amplitude is for a dipole of transverse size x scattering on the target nucleus at 
the impact parameter b. p is the density of nucleons in the nucleus, T{b) is the nuclear profile 
function equal to the length of the nuclear medium at the impact parameter b, such that 
T(6) = — IP' for a spherical nucleus of radius R. Also, x±_ = \x\, Cp = ^2N^ ^ 

the (non-perturbative) momentum scale characterizing each nucleon. We have put as the 
bare coupling constant to underline that the calculation of [65] leading to Eq. flAl|) was done 
for fixed coupling. Until the running coupling corrections are included the coupling in Eq. flAll) 
is the bare coupling. 

One can show (see [30] for a pedagogical derivation) that the expression in the exponent of 
Eq. ( lAll) is equal to 

_pr(6)^ (A2) 



with 0"''''^ the scattering cross section of the dipole on a single nucleon. The dipole-nucleon 
cross section due to a two-gluon exchange is given by 



= 2 «J ^ / (2 - - e--) = 2 «J ^ vr .Mn -i-, (A3) 



^ j [Z2]2 V ' \x\k 
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where the two factors of in the denominator come from the propagators of the two gluons. 
A is used as the infrared cutoff of the /-integral in Eq. (lA3p . 

To include running coupling corrections into Eq. (lAip . it is sufficient to include them into 
the cross section (IA3p . The corresponding "dressed" diagram is shown in Fig. [161 where the 




Figure 16: The interaction of a qq dipole with a single nucleon in the nucleus in the Glauber- 
Mueller approximation "dressed" with quark bubble corrections. 

two exchanged gluons have chains of quark bubbles on them. Indeed one should sum over 
connections of each of the gluons to the quark and the anti-quark in the dipole: this is actually 
what gives the factor (2 — e*'''^ — e~*'''") in Eq. (lA3p . Since in the eikonal approximation the 
gluons light cone momentum is small, /_ = 0, the gluon virtuality is = —t^. The quark 
bubble chains introduce two usual denominators into the fixed coupling equation (lASp . turning 
it into 

2 (2 - e*'-^ - e^*'-^) , (A4) 

where we have completed the factors of Nj to the full QCD one-loop beta-function, Nf 
—6 TT (32- Eq. flA4p obviously contains the physical running coupling squared and can be rewrit- 
ten as 

a--^ = 2^1^ al{f) (2 - e^'- - e^'-) . (A5) 

Fourier-transform of a running coupling constant, as is implied in Eq. (lASp . is always dan- 
gerous due to Landau singularity. Avoiding uncertainties introduced by Landau pole we assume 
that I a; I ^ 1/A ^ 1/Aqcd and approximate the integral in Eq. flASP by 

i/|x| 

A 

where we have made use of Eq. fl46l) . While indeed in real life A ~ Aqcd j we have made the A ^ 
Aqcd assumption to be able to disentangle perturbative and non-perturbative contributions in 
the integral of Eq. flASp . One should then take the A — »• Aqcd limit to identify the non- 
perturbative part of the expression. An analogy would be a finite temperature medium, where 



qqN 



77212 



1 + /^s In ^ 
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A ~ T and for high temperatures T 3> Aqcd our above approximation would be justified. 
We would then obtain a factor of ^^(T^) in Eq. (lA6p . which would become a non-perturbative 
factor for lower temperatures T ~ Aqcd- 

Comparing Eq. flA6l) with Eq. flASp we see that the inclusion of running coupling corrections 
is accomplished by the replacement 

"J (^) (^') • (A"^) 

Eq. dSD with the running coupling corrections included becomes 



iV(x, h) = I- exp 



For each multiple rescattering one of the coupling constants is determined by the non- 
perturbative scale A and could be large. That is why this coupling constant is sometimes 
absorbed into the gluon distribution function of the nucleon [48,65,81-83]. The other coupling 
comes with the scale l/|a;|, which, for perturbatively small dipoles (corresponding to production 
of particles with large transverse momenta) is large, making the corresponding coupling constant 
small. 
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